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Abstract. Long-lasting, intense, stellar X-ray flares may 
approach conditions of breaking magnetic confinement 
and evolving in open space. In the perspective of search¬ 
ing for possible tracers of non-confinement, we explore this 
hypothesis with hydrodynamic simulations of flares occur¬ 
ring in a non-confined corona: model flares are triggered by 
a transient impulsive heating injected in a plane-parallel 
stratified corona. The plasma evolution is described by 
means of a numerical 2-D model in cylindrical geometry 
R, Z. We explore the space of fundamental parameters. As 
a reference model, we consider a flare triggered by a heat¬ 
ing pulse of 10 erg cm“^ s“^ lasting 150 s and released 
in a region ^ 10® cm wide and at a height ^ 2 x 10® 
cm from the base of the stellar surface. The pressure at 
the base of the corona of the unperturbed atmosphere is 
0.1 dyne cm“®. The heating would cause a 20 MK flare 
if delivered in a 40000 km long closed loop. The modeled 
plasma evolution in the heating phase involves the prop¬ 
agation of a 10 MK conduction front and the evaporation 
of a shocked bow density front upwards from the chro¬ 
mosphere. As the heating is switched off, the temperature 
drops in few seconds while the density front still prop¬ 
agates, expanding, and gradually weakening. This kind 
of evolution is shared by other simulations with different 
coronal initial pressure, and location, duration and inten¬ 
sity of the heating. The X-ray emission, spectra and light 
curves at the ASCA/SIS focal plan, and in two intense 
X-ray lines (Mg XI at 9.169 A and Fe XXI at 128.752 A), 
have been synthesized from the models. The results are 
discussed and compared to features of confined events, 
and scaling laws are derived. The light curves invariably 
show a very rapid rise, a constant phase as long as the con¬ 
stant heating is on, and then a very fast decay, on time 
scales of few seconds, followed by a more gradual one (few 
minutes). We show that this evolution of the emission, 
and especially the fast decay, together with other poten¬ 
tially observable effects, are intrinsic to the assumption of 
non-confinement. Their lack indicates that observed long- 
lasting stellar X-ray flares should involve plasma strongly 
confined by magnetic fields. 
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1. Introduction 

Stellar coronae are not spatially resolved with present day 
observations. It is widely accepted that they share many of 
the basic mechanisms of the solar corona, the only resolved 
corona that we know. Decades of X-ray observations have 
shown that both its structuring and heating are dominated 
by the ambient magnetic field. Evidence has been found, 
such as the presence of coronal plasma at temperatures 
of several million degrees (e.g. Linsky 1981, Schmitt et 
al. 1990), that some form of confinement of X-ray bright 
plasma should be at work on stars, probably associated 
with intense magnetic fields. It is then natural to wonder 
what is the morphology of the confining structures in a 
stellar corona. On the Sun, plasma is typically confined in 
magnetic arches (the loops), and both the geometry and 
luminosity of many of them are steady for relatively long 
times. 

If stellar coronae are similar to the solar one, then 
structures similar to solar loops should be found on stars. 
Several approaches are possible to determine the dimen¬ 
sions of coronal structures. An eclipse on Algol has al¬ 
lowed to find the volume responsible of a flare (Schmitt 
& Favata 1999). Loop fitting has been used to constrain 
the average linear dimensions of stellar coronal loops (e.g. 
Maggio & Peres 1997). Another source of information on 
the characteristic of stellar coronal loops is the evolution 
of X-ray flares. In particular, basic mechanisms of plasma 
cooling are radiation and thermal conduction from the 
corona downwards to the cooler chromosphere. Linking 
the plasma characteristic cooling times to the observed 
flare decay time allows to infer the linear size of the cool¬ 
ing region. This property has been largely exploited in 
the past to estimate the size of stellar flaring loops (e.g. 
Haisch 1983, van den Oord et al. 1988, van den Oord & 
Mewe 1989, Reale et al. 1988, Reale et al. 1997, Favata & 
Schmitt 1999, see Reale 2001 for a review). 

One implication of the simple scaling from cooling 
times is that slower and slower decays, associated with 
long-lasting flares, correspond to longer and longer loops. 
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As a matter of fact, stellar X-ray flares are typically ob¬ 
served to last much longer than average solar X-ray flares, 
up to several days whereas the duration of solar flares 
ranges from few minutes to several hours. In spite of the 
obvious observational bias toward detecting long events, 
the evidence of long flares has been taken as strong indica¬ 
tion of very large flaring regions, in particular comparable 
to, or even larger than, the stellar radius (e.g. Graffagnino 
et al. 1995, Favata & Schmitt 1999, Tsuboi et al. 2000, 
Favata 2001). On the Sun, large scale and long-duration 
events typically occur in complex coronal structures, in 
which the magnetic field undergoes major rearrangements 
and several loop structures are progressively involved (e.g. 
two-ribbon flares): the large amounts of energy released 
lead to large plasma pressure which may exceed the mag¬ 
netic one and cause the end of confinement. The ques¬ 
tion now is: may plasma in long-lasting and intense stellar 
flares break the magnetic cage, and erupting, be free to 
move and evolve in open space? 

This work addresses this question, by modeling flares 
occurring in a totally non-conflned stratified corona. Al¬ 
though we do not pretend this to be an entirely realistic 
scenario, it can be considered as an extreme limit of situ¬ 
ations in which the plasma dynamics dominates over the 
confining magnetic field and therefore governs the X-ray 
luminosity evolution. 

Indeed, some solar and stellar flares show phenomena 
that indicate that at least some of the plasma is not con¬ 
fined. In the solar case, coronal mass ejections (CMEs) 
following flares are not confined. In the stellar case VLBI 
measurements of nonthermal particles propagating away 
from the star are observed. A good example is a flare in UX 
Ari observed by Mutel et al (1985). White and Franciosini 
(1995) also show that the relativistic plasma responsible 
for the gyrosynchrotron radio emission in stellar comes 
from expanding plasma. 

Here we focus on studying the X-ray emitting plasma. 
Our scope is to explore: to what point actual stellar flares 
can approximate free-expanding flaring plasma, if there 
are distinctive signatures in X-ray emission (obtained from 
the modeling) which characterize unconfined flares vs con- 
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fined ones and therefore allow us to discriminate the con¬ 
fining role of the magnetic field, and if observed stellar 
flares present such signatures. 

The study that we illustrate here presents also interest¬ 
ing results on the theoretical point of view: we describe the 
evolution of the non-confined plasma in a stratified coro¬ 
nal atmosphere, i.e. which is hot ab initio and where ther¬ 
mal conduction is fully efficient. This aspect makes this 
work quite different from typical models of bursts, like su¬ 
pernovae, which instead propagate in non-conducting and 
cooler media. 

Our approach is to set up an initial stratified atmo¬ 
sphere, similar to a confined one, including chromosphere, 
transition region and corona, and where the plasma is not 
confined to move along one spatial direction, and then to 
impart a localized energy impulse. Plasma dynamics and 
heat conduction can occur isotropically. The evolution of 
such system needs to be described by a numerical 2-D 
hydrodynamic model, with enhanced spatial resolution in 
the low corona and in the transition region and including 
plasma thermal conduction; therefore it is also a demand¬ 
ing numerical problem. In order to compare our results 
with observations and to obtain diagnostical feedbacks, we 
synthesize from our model results the expected evolution 
of the integrated X-ray flare luminosity, i.e. by simulat¬ 
ing a flaring unresolved source, such as a star. To have a 
realistic scenario, representative of real flare observations, 
we have chosen to synthesize the emission in a wide X-ray 
band such as at the focal plane of the ASCA/SIS (sen¬ 
sitive in the range, approximately, 0.3 - 10 keV), and in 
two intense X-ray lines. Similar results can be expected for 
other wide band CCD instruments such as Chandra/ACIS 
and XMM/EPIC, while the lines may be detectable by in¬ 
struments with good spectral resolution such as Chandra 
and XMM-Newton grating detectors. 

In Sect. 0 we describe our modeling approach in de¬ 
tail, in Sect.^ the simulations performed and the results 
obtained are presented, in Sect. ^ the results are discussed 
and in Sect. conclusions are drawn. 

2. Modeling 

2.1. The Model Set Up 

The general concept of our modeling non-confined coronal 
flares is similar to modeling confined flares (e.g. Peres et 
al. 1982): the flare is triggered by an impulsive heating 
injected in a stratified plane-parallel static corona. The 
basic difference is that the unperturbed atmosphere is not 
a semicircular coronal loop and that the heat transport, 
as well as any plasma motion, is allowed in any direction, 
instead of being limited only along magnetic field lines 
parallel to the loop flux tube. The evolution of a heat 
perturbation propagating in any direction in a stratified 
atmosphere can be properly described with a 2-D cylindri¬ 
cal geometry {R,Z). This description is adequate for our 


purposes of deriving the global properties of the plasma 
evolution and X-ray emission as it would be detected at 
stellar distances. 

The solar template of confined coronal flares has 
guided us in configuring the initial unperturbed atmo¬ 
sphere and in tuning the plasma and heating parameters. 
As shown in Fig. the initial atmosphere consists of a 
stratified plane-parallel corona extending vertically from 
the stellar surface. The limit of our computational do¬ 
main is at 10^^ cm from the stellar surface. The corona 
is maintained steady at coronal conditions by a uniformly 
distributed heating. As our reference starting conditions 
we have taken an atmosphere with a pressure at the base 
of the corona of 0.1 dyne cm“^, typical of coronal quiet 
regions, which may reasonably approximate non-confined 
coronal conditions. In these conditions the temperature of 
the upper layers of the non-flaring corona of our compu¬ 
tational domain is below 3 MK. We have made some sim¬ 
ulations also with a base pressure 1 dyne cm“^, a denser 
and hotter (< 6 MK) initial atmosphere, to check the ef¬ 
fect of changing the initial conditions to rather extreme 
cases of very active ” non-confined” coronae. The bottom 
of the corona is matched to a much cooler chromosphere 
through a steep transition region. All this setting is similar 
to that for modeling confined coronal flares (e.g. Peres et 
al. 1982, Cheng et al. 1983, McClymont & Canfield 1983). 
The chromosphere is required as mass reservoir for the 
flare evolution; our choice has been to consider an isother¬ 
mal (2 X 10^ K) and non-radiating chromosphere (Cheng 
et al. 1983). While the detailed structure of the chromo¬ 
sphere is important for the fine details of the evaporation 
and coronal evolution, a coarse chromospheric model like 
the one used here suffices to show the general characteris¬ 
tics of the evolution addressed here. We have verified that 
the unperturbed atmosphere is steady and stable for times 
much longer than the flare evolution times considered for 
this study. 

The plasma evolution is described by solving the 2- 
D plasma hydrodynamic equations of conservation of 
mass, momentum and energy with the geometry described 
above: 

^-^V-(nv)=0, (1) 

^ + V ■ (nv V ) = Vp -H ng , (2) 

at m 

du 

— -|-V-[(u -|-p)v ] = mng v -n'^P{T)+'H{r )-VFc , (3) 

u=lmnv^+£, p = , T=^, F, = 

2 2kn 

where n particle density per unit volume, t time, v the 
plasma bulk velocity, p pressure, m average mass per par¬ 
ticle (to = 2.1 10“^'^ g for solar abundance), g gravity 
acceleration, P(T) radiative losses per unit emission mea¬ 
sure, TL heating per unit volume (uniform and constant 
in transition region), £ is the internal energy density, T 
plasma temperature, heat flux, k plasma thermal con- 
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Fig. 1. Initial configuration of the reference model (cf. model 1 in Table 1); the geometry is cylindrical with coordinates 
r,z. Panels (a) and (b) show the temperature and density distributions along the vertical z axis (r = 0). Panel (c) 
shows the spatial distribution of the impulsive heating vs radial distance (R) and height (Z), where the radius of the 
dark grey circle is an (see Eqs. Si- Notice that the upward vertical direction is rightwards in panel (c). Panels (a), 
(b) and (c) show a limited region of the whole computational domain. Panel (d) shows the temperature and density 
(dashed) distributions along the vertical z direction on a much larger scale. 
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ductivity and 7 = 5/3. The flare is assumed to occur on a 
solar-like star, with solar surface gravity and radius. 

Plasma thermal conduction is isotropic. Radiative 
losses are those of Raymond and Smith (1977) but are set 
to zero for T « 2 X 10^'’ K, i.e. for chromospheric plasma. 
The heating term generally consists of a steady heat¬ 
ing term which maintains the unperturbed atmosphere in 
thermal equilibrium and of a transient term, which trig¬ 
gers the flare. In our simulations, the steady heating has 
been set to zero both in the chromosphere and elsewhere, 
since the transient heating is switched on, in order to de¬ 
scribe the free plasma cooling with no interference from 
any heating source. 

The transient heating term has been assumed as a sep¬ 
arable function of space and time: 

H{Z,R,t)^Hof{Z,R)g{t) (5) 

The spatial distribution of the heating is a 2-D circular 
Gaussian: 


f{Z,R) = exp 


[Z- Znf + {R-RHf 




( 6 ) 


As for the temporal evolution g(t) the impulsive heat¬ 
ing has been assumed to be switched on at t =0 (g( 0 ) = 0 ), 
kept constant {g{t) = 1 ) for a given time lapse 0 < t < Ih 
and then switched oS g{t > tn) = 0- 


2.2. The Code 

The conservation equations shown above are solved nu¬ 
merically on a fixed geometrical grid. The time integration 
is made with “time splitting” in two separate phases: a) 
the proper hydrodynamic phase, which includes the mass 
and momentum equations and the advective terms of the 
energy equation; b) the thermal conduction phase, which 
includes also radiation and heating processes. The hydro- 
dynamic section is based on an explicit Flux Corrected 
Transport (FCT) numerical scheme, improved for effi¬ 
ciency and accuracy in large-scale smooth regions (Reale 
et al. 1990). The highly non-linear thermal conduction 
is treated with an Alternating Direction Implicit (ADI) 
scheme (Reale 1995), which warrants unconditional sta¬ 
bility. The two phases are solved independently according 
to a time-splitting approach, i.e. using different integra¬ 
tion time steps, which depend on the characteristic times 
of the physical effects involved in each section and on the 
stability criteria of the respective numerical schemes, i.e. 
FCT and ADI. 

As mentioned above the geometrical domain is 2-D in 
cylindrical coordinates {R,Z). The grid does not change 
during a simulation and the grid cells are rectangular: the 
spacing dZ along the Z direction is uniform in the chromo¬ 
sphere {dZ = 2x 10® cm) and then monotonically increases 
with a geometric progression up to dZ ~ 2 x 10 ® cm at 
Z « 10^^ cm; this spacing warrants a stable and steady 
atmosphere in equilibrium conditions; the spacing along 
R is uniform as far as 7 x 10^® cm from the central axis 
and then monotonically increases as far as ^ 10^^ cm. The 


spacing along Z and R has been chosen so as to have a 
relatively high resolution and a regular grid in the internal 
domain region, where the energy is released and the flare 
mostly occurs, and at the same time, a very large domain 
(Bedogni & Woodward 1990), so that the boundary con¬ 
ditions at the far extremes become unimportant for the 
problem under analysis. The space cells along Z are 800 
for all simulations performed. Along R most simulations 
has been done with 200 grid cells for the internal uniform 
region, plus 28 cells of non-uniform expanding grid; some 
simulations have been done with 100 and 400 grid cells for 
the internal region. 

Reflecting boundary conditions have been set for the 
central Z (R = 0) axis. Fixed density and temperature 
have been set constant on the lower R axis {Z = 0), and 
with zero gradient on the upper one. The Z component 
of the velocity is set to zero at both Z boundaries, the R 
component with zero gradient. At the far extreme of the 
R domain zero gradient conditions have been set for all 
the relevant quantities. 

The code has been parallelized on the basis of a geo¬ 
metric domain decomposition and using the High Perfor¬ 
mance Fortran (HPF). In the Alternating Direction Im¬ 
plicit (ADI) scheme, the equation is integrated at each 
time step in two half steps. At each half-step the numer¬ 
ical integration advances in one of the two spatial direc¬ 
tions. The integration scheme at each half-step is implicit: 
the solution at one grid-point is derived from the solution 
at the same time step at an adjacent point. The data de¬ 
pendencies prevent direct parallelization of the nested DO 
loops, as in the case of the FCT section: we have adopted 
a parallelization solution which minimizes the interproces¬ 
sor communications and the sequential operations, with¬ 
out changing the numerical algorithm. Simulations have 
been executed on the Cray T3E at CINECA (Bologna, 
Italy), on a Compaq server and on a Compaq HPC Clus¬ 
ter at Osservatorio Astronomico di Palermo (Italy). 


3. The Simulations 


3.1. The Parameter Space and The Reference Flare 

Table ^ shows the parameters and some results of the sim¬ 
ulations of non-confined flares presented here (plus one 
model of confined flare for comparison). For each simu¬ 
lation, there are the basic simulation parameters, i.e. the 
pressure at the base of the initial unperturbed corona, and 
the parameters of the heating release (see Eqs. || and ||), 
i.e. the heating rate per unit volume Hq, the height Zh of 
the heating maximum, the Gaussian width an, the total 
rate, and the duration of the heating (Eq. ^), and some 
relevant results, i.e. the plasma maximum temperature, 
the maximum density at the position of the maximum 
temperature, the maximum velocity throughout the com¬ 
puted flare evolution, and, finally, the light curve deca y 
time after the initial very fast decay phase (see Sect. 3.3). 
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Our simulation strategy is to set up and compute a 
reference model case and then perform other calculations 
by changing the parameters relevant for our study one 
by one. Our reference configuration has a coronal base 
pressure po = 0.1 dyne cm“^, which is a value typical of 
quiet coronal regions or even of open coronal structures 
on the Sun. Since we suppose that higher pressures may 
occur on active stars, values up to po = 1 dyne cm“^ have 
also been considered. 

In our reference simulation the flare is triggered by 
a heating pulse with maximum rate rate per unit vol¬ 
ume Hq = 10 erg cm“^ s“^ and width of the Gaussian 
an = 5 X 10® cm centered on the Z axis {Rh = 0) at a 
height oi Zh = 2x 10^ cm above the base of the atmo¬ 
sphere, well above the transition region (Eq. Similar 
heating rate, width and position would, for instance, ig¬ 
nite a 20 MK confined solar flare at the top of an active 
region loop (Peres et al. 1987), which may be taken for 
comparison to the results presented here. As shown in Ta¬ 
ble 1^, other simulations have been performed with higher 
values of the rate Hq up to 100 erg cm“® s“^, in order 
to explore the trend toward more intense events such as 
those typical of active stars. We have also considered heat 
deposited higher in the corona, i.e. more distant from the 
stellar surface, and, at the same time, over a more ex¬ 
tended volume, while either maintaining the same total 
heating rate or doubling it. In order to obtain a signifi¬ 
cant brightening over the background atmosphere, in all 
simulations with higher base coronal pressure, the heating 
rate Hq is higher than the reference value, since the atmo¬ 
sphere is brighter and also more efficient in radiating the 
additional heat, this making any smaller heating barely 
effective. 

For most simulations, the duration of the constant 
heating (the time in which g{t) = I) is 150 s, but a case of 
much longer heating {tn = 600 s) has also been explored. 

The reference simulation (named re/in Table will be 
described in more detail in the following. With the total 
heating rate of ~ 2 x 10^® erg s“^, the total heat released 
is therefore ~ 3 x 10®® erg. The evolution is computed for 
a total time of 800 s, i.e. the decay is studied for 650 s. 
The grid is 800 x 228 grid cells. 

For comparison with the models of flares in non- 
confined atmospheres. Table |l| shows also the parameters 
and results of a model of flare in a confined coronal loop, 
labelled with Conf. In particular, the parameters are those 
of a flare in an active region loop of half-length 2 x 10® 
cm and initial pressure 6 dyne cm“® (Peres et al. 1982, 
Peres et al. 1987). The heating duration, spatial extent 
and intensity per unit volume have been chosen identical 
to those of the reference model. A loop cross-section area 
of 2.5 X 10^^ cm® has been assumed, corresponding to a 
loop aspect, i.e. radius divided by half-length, of 0.1. 

3.2. The plasma evolution of the reference flare 


3.2.1. The heating phase 

Fig. ^ shows distributions of temperature, density, pres¬ 
sure and vertical component of velocity along the central 
Z axis at various times during the phase in which the 
heating is switched on in the reference simulation. Fig. 
shows the analogous results of the evolution after the heat¬ 
ing has been switched off. Figs. and || respectively show 
grey-scale images of the temperature and density contrasts 
T/Tq and n/ng, where Tq and uq are the temperature and 
density of the unperturbed atmosphere at various times 
during the same simulation, including both heating and 
decay phase. Figs. ^ and ^ clearly show the evolution 
of the plasma in the heating phase. 

Because of the strong impulsive heating the plasma 
temperature at the center of the heating release on the 
Z axis increases rapidly (in less than 2 s) from 1.5 to 12 
MK (and similarly the pressure by a factor 10) and then 
remains practically constant there. A thermal conduction 
front develops and propagates spherically and very rapidly 
in such a thin corona, taking ~ 2 s to reach the base of 
the corona, the transition region and the chromosphere, 
and, on the opposite side, a height Z ~ 5 x 10® cm. As the 
conduction front propagates in corona, expanding from 
the heating region, it fades and slows down significantly. 
While a height Z = 10®® cm is reached by the temperature 
change after only ^ 15 s, it takes one minute more to reach 
Z = 2x 10®® cm, and at the end of the heating the conduc¬ 
tion front is only slightly higher than that. In the orthog¬ 
onal (radial) direction the conduction front is slower, be¬ 
cause the temperature increases upwards, and therefore so 
does the efficiency of thermal conduction, while this does 
not occur in the radial direction; R = 10®® cm is reached 
only after 50 s since the beginning of the simulation and 
at the end of the heating the distance of the conduction 
front from the Z axis is i? = 1.5 x 10®® cm. To have an 
idea of the weakening of the expanding conduction front, 
note that, at the end of the heating, the temperature is 
below 10 MK farther than 2 x 10® cm from the center of 
the heating release and below 4 MK for R > 10®® cm. 

As soon as the chromosphere is hit and heated by the 
conduction front, it expands upwards with a shocked den¬ 
sity front^. This is the chromospheric evaporation modeled 
in confined solar coronal flares (e.g. Nagai 1980, Peres et 
al. 1982). The plasma maximum speed of the evaporation 
front rapidly increases to ~ 600 km/s in about 30 s and 
then remains between 600 and 700 km/s. The pressure of 
the evaporation front along the Z axis increases by an¬ 
other factor 5, reaching a value ~ 10 dyne cm“® at the 
base of the corona. The density front reaches the central 
position of the heating along the Z axis in ^ 30 s, and a 
height Z ^ 10®® cm at the end of the heating (t = 150 s), 
corresponding to an average propagation speed of ^ 600 

® The code does not include the effect of plasma viscosity, 
which has been shown to be important in solar flares and may 
spread out and therefore inhibit the shock front. 
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Table 1. Parameters of simulated non-confined flares -I- one confined flare 


Name 

pS 

dyn/cm^ 

Hg 

erg cm ® s ^ 

10® cm 

d 

10® cm 

mot 

10®®erg/s 

S 

'T'fl 

max 

10’’ K 

10® cm"® 

* max 

km/s 

' dec 

S 

Ref 

0.1 

10 

2.0 

0.5 

2 

150 

1.3 

0.9 

800 

290 

Long 

0.1 

10 

2.0 

0.5 

2 

600 

1.3 

0.9 

1000 

380 

H30 

0.1 

30 

2.0 

0.5 

6 

150 

1.7 

1.2 

1000 

270 

HlOO 

0.1 

100 

2.0 

0.5 

20 

150 

2.5 

1.9 

1200 

320 

X5H5 

0.1 

5 

5.0 

1.0 

8 

150 

1.6 

0.3 

700 

180 

X5H10 

0.1 

10 

5.0 

1.0 

16 

150 

2.0 

0.4 

900 

210 

P03H30 

0.3 

30 

2.0 

0.5 

6 

150 

1.7 

1.7 

1000 

350 

P03H100 

0.3 

100 

2.0 

0.5 

20 

150 

2.5 

2.6 

800 

370 

PI 

1.0 

100 

2.0 

0.5 

20 

150 

1.9 

2.8 

600 

200 

“Con? 

6 

10 

2.0 

0.5 

0.2 

150 

2.1 

75 

360 

250' 


“ Pressure at the base of the corona 
** Maximum heating rate per unit volume (see Eqi 
Height of the center of the heating (see Eq-i) 

Width of the heating (see Eq|) 

® Total heating rate 
^ Heating duration (see Eq.|) 

® Flare maximum temperature 

^ Maximum density at the position of the temperature maximum 

* Maximum velocity 

^ Light curve (in the ASCA/SIS band) e-folding decay time after the initial fast decay phase 
^ Model of a flare confined in a solar active-region loop 

* E-folding decay time since the light curve maximum 


km/s. As a consequence the density in the central post¬ 
shock region increases above a factor 3 but practically in 
all the coronal region it is < 10® cm“^ at any time. 

Given the highly localized heating release in a non 
confined atmosphere, the evaporation front, which mostly 
determines the soft X-ray burst, is not plane-parallel: it 
is practically semicircular, stronger along the Z axis and 
weaker and weaker moving farther and farther from R = 0. 
The reason is that the conduction front arrives at the chro¬ 
mosphere first along the central Z axis, and progressively 
later at greater and greater distances from it. The front 
amplitude gets weaker and weaker and with a smaller 
pitch angle. For R > 10^® cm the evaporation speed is 
< 100 km/s and the local density enhancement very small 
at any time. 

The density in the post-shock region is not uniform: 
the plasma behind the shock moves faster than at the 
shock and accumulates there; the evaporation front be¬ 
comes then a bow front and a region of low density (a 
proper depression) forms behind it. The bow front fades 
farther and farther from the central head. The diameter of 
the front is less than 2 x 10^® cm during the whole heating 
phase. 

3.2.2. The decay phase 

Fig.^ shows distributions of temperature, density, pressure 
and vertical component of velocity along the central Z 


axis at various times of the decay phase, in the reference 
simulation. 

As the heating is switched off (t = 150 s), the temper¬ 
ature suddenly decreases in the heating region: it halves 
(from 12 to 6 MK) in about 5 s at the center of the heating 
release. The cause is the cooling by thermal conduction, 
particularly efficient because the density of the hot plasma 
remains relatively low (< 10® cm“®). The cooling front 
propagates radially from the heating region, while the in¬ 
ternal temperature gets lower and lower. The hot front 
produced by the heating impulse has instead weakened 
and practically disappeared when reaching the higher and 
relatively hotter (but very thin) corona. After one minute 
since the heating switch off (t « 210 s), the temperature 
along the Z axis overshoots below the initial temperature 
of ~ 1 MK. At t = 300 s only a shell at ^ 10^® cm from 
the system origin, and ^ 10^® cm thick, is still (slightly) 
hotter than the initial atmosphere, while the originally 
heated region is practically all cooler than 1 MK. 

While such a fast cooling is occurring in the heated 
region, the evaporation front continues to propagate up¬ 
wards and outwards. The shell expands but maintains 
more or less the same shape, becoming geometrically 
thicker and thicker. The head front is at speed around 
400 km/s, the back front is steadily slightly below 800 
km/s. The difference in speed causes the accumulation. 

The depression formed behind the head front deepens 
and expands continuously behind the evaporation front 
after the heating is switched off. At first it is confined to a 
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Z(10® cm) 


Fig. 2. Distributions of temperature T, density n, pressure p, and vertical component of velocity V = Vz, along the 
central Z axis (i? = 0) and at various times (s) during the heating phase of the reference simulation of the non-confined 
flare. The arrows indicate the position of the center of the heating distribution. 


rather small region behind the head front, but it extends 
for ~ 10^° cm below the head front at t = 300 s. 

At this time, the thickness of the head front, with a 
density ~ 5 times higher than the background atmosphere, 
is around 5 x 10® cm; the front has reached a height Z ~ 
2 X 10^® cm above the stellar surface and has a diameter 
(distance from the R = 0 axis) of ~ 1.5 x 10^® cm. The 
density in the head front reduces from ~ 3 x 10® cm“® 
to ~ 10® cm“® as it moves away, while the region behind 
gets less and less dense, the core of it decreasing below 
10^ cm“®. The pressure in that core, but even below it, 
in a layer 10^® cm thick, is around 0.01 dyne cm“®, 1/10 
smaller than the initial one. 

After t = 300 s the cooling of the central regions and 
the expansion of the evaporation shell progress with no 
new feature to remark. At t = 800 s no temperature en¬ 
hancement is visible anywhere, while the evaporation front 
is still well evident with the top above 3 x 10^® cm from 
the surface, a height where the density of the background 


corona is below 10^ cm“®, and the tail all above 10^® cm. 
The diameter of the shell is ~ 5 x 10^® cm, its maximum 
thickness ~ 1.5 x 10^® cm, but the density is not uniform 
inside, ranging between 3 x 10^ cm“® (ahead) and 10^ 
cm“® (behind). The speed of the head front is decreas¬ 
ing to 200 km/s. The gap behind the front appears to be 
gradually Ailing and reducing. 

3.3. The X-ray emission of the reference flare 

From the model results we have synthesized the X-ray 
emission both in the ASCA/SIS band and in two rele¬ 
vant X-ray lines, namely Mg XI at 9.169 A and Fe XXI at 
128.752 A with peak formation temperatures of 6.5 MK 
and 10.5 MK, respectively. These two lines are represen¬ 
tative of strong emission lines typically visible in the band 
of grating detectors such as those on board Chandra and 
XMM-Newton and with formation temperatures relevant 
to become significant lines during the reference model. 
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Fig. 3. As in Fi^ but taken at times after the heating is switched off. 


In the longer wavelength line the spectrometers typically 
have higher spectral resolution. 

The X-ray spectra at the focal plane of the ASCA/SIS 
and in the two selected lines, have been computed in the 
whole computational domain and at each output time, 
according to: 

£es{E, r, z, t) = n^{r, z, t) G{E, T{r, z, t)) (7) 

where G{E, T) is the emissivity per unit emission measure 
and per unit photon energy in the case of the ASCA/SIS, 
and per unit Angstrom in the case of the spectral lines. 
The emissivity in the ASCA/SIS band from 0.6 to 14 
keV has been obtained by folding the ASCA/SIS effective 
area and response with Raymond & Smith (1977) spectra 
of isothermal plasma, assuming full Anders & Grevesse 
metal abundances and a small interstellar medium ab¬ 
sorption {Nh = 10^® cm“^). The emissivity in the two 
spectral lines have been extracted from the CHIANTI 
database (version 3.02, Dere et al. 1997), assuming ion¬ 
ization equilibrium computed according to Arnaud and 
Raymond (1992) and solar metal abundances (Grevesse 


and Anders 1991) and using the default density of 10^° 
cm“^. 

The spatial distribution of the total focal-plane emis¬ 
sion per unit volume at each time is then obtained by 
integrating in energy: 

£s{r,z,t)=f £es{E,r,z,t) dE . (8) 

^ Emin 

Fig. I shows, for the reference simulation, cross- 
sections of the resulting distributions of plasma emission 
in the ASGA/SIS band at the same times and in the same 
region as Figs.|| and ||. The cross-sections are taken on 
a plane across the Z axis. The emission at early times is 
clearly concentrated in the heating region and mostly in 
a thin layer at the base of corona. Until t = 30 s an al¬ 
most spherical weak emission front expands progressively 
and is clearly associated to the propagation of the ther¬ 
mal front (see Fig.^. After t = 30 s another brighter front 
rises from the base of the corona, with a less regular shape, 
and is due to the density front evaporating from the chro¬ 
mosphere (see Fig.||). Soon after the heating is switched 
off (t = 154 s), the bright knot at the base of the corona 
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Fig. 4. Grey scale images of the temperature contrast T/Tq, where To is the temperature distribution of the unper¬ 
turbed atmosphere, and velocity held {arrows) at various times during the reference simulation. It is shown the region 
around the heating deposition. The grey-scale is logarithmic. The maximum speed corresponding to the longest arrow 
is shown below each image. 


rapidly fades, d ue to the sudden temperature drop de¬ 
scribed in Sect. 3.2.2 . At t = 200 s most of the emission 
comes from the evaporating front, which continues to ex¬ 
pands but gradually fades away (t = 300 s). 


Since the hare would not be resolved at stellar dis¬ 
tances, we have then integrated the emission distribution 
over the relevant region of the computational domain to 


obtain the light curve at the focal plane of ASCA/SIS and 
the hux at Earth in the two selected spectral lines: 

£{t) = [ £s{r,z,t) dV (9) 

Jv 

where V — irR^^^Z^ax is the domain volume. In deriv¬ 
ing the focal-plane light curve £{t) we have assumed a 
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Fig. 5. As in Fig.|| but for the density contrast njno, where no is the density distribution of the unperturbed atmo¬ 
sphere. 


standard distance of the flaring source of 1 pc from the 
detector. 

Figure ^ shows the flare net light curves at the 
ASCA/SIS focal plane, obtained by subtracting the count 
rate at the initial time of each simulation, taken as back¬ 
ground non-flaring emission, for all cases reported in Ta¬ 
ble 1^, including the model of confined flare. 

Let’s first focus the attention to the light curve of the 
reference case (the lowest one) and compare it with the 


confined one: it looks quite different. The light curve of 
the confined flare is characterized by a relatively rapid 
rise phase, a well-defined peak and a gradual decay. The 
light curve of the modeled unconfined flare, instead, after a 
very rapid rise (~ 20 s) to the maximum value 5 times 
the initial emission level), shows a long plateau, staying 
virtually constant, with small fluctuations 20%), until 
the end of the heating (150 s). As soon as the heating is 
switched off, the light curve suddenly drops: the emission 
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Fig. 6. Cross-sections (grey scale) of the X-ray emission (logarithmic scale normalized to global maximum) distribution 
in the ASCA/SIS band on a plane across the Z-axis and in the region where the heating is released, and velocity field 
(arrows) at various times during the reference simulation. Region, grey scale and maximum speed as in Fig^. 


halves in less than one minute. We will call this phase 
Fast decay. Since then on, the emission decreases more 
gradually with an e-folding decay time of ^ 5 min, as 
reported in Table |^. We will call this phase Gradual decay. 
We notice that the light curve of the non-confined flare 
follows much more closely the evolution of the heating: in 
particular, the emission drops as soon as the heating stops 
whereas, in the confined case, it keeps on increasing and 


begins to decay only ~ 200 s after the end of the heating. 
The decay trends are complementary: in the confined case 
it is slow at first and then gradually becomes faster and 
faster, also because more and more plasma cools and is no 
longer visible in the instrument passband. 

We have made the exercise of a fit of ASCA/SIS 
spectra sampled along the light curve with isothermal 
plasma spectra, a standard analysis in stellar coronal 
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Time (s) 


Fig. 7. Light curves at the ASCA/SIS focal plane, assuming a distance of 1 pc, synthesized from the results of the 
simulations listed in Table The light curve of a modelled confined flare is also shown for comparison {thickest solid 
line, labelled CONF). 


physics. We obtain a sequence of temperature and nor¬ 
malization (emission measure) values which can be put in 
a density-temperature diagram (where the square root of 
the emission measure is used as proxy of the density, as 
it has been typically done for confined flares[|, e.g. Syl- 
wester et al. 1993, Reale et al. 1997), in order to compare 
this work with previous work to model stellar flare decay 
(Reale et al. 1997, Reale & Micela 1998, Favata & Schmitt 
1999). This diagram is used as a diagnostics of the possible 
presence of sustained heating during the decay (Sylwester 
et al. 1993). 

The upper panel of Fig. || shows the path of the ref¬ 
erence non-confined flare in the density-temperature dia- 


^ Notice that for non-confined flaring plasma, this assump¬ 
tion, based on constant volume, is grossly incorrect. 


gram, and, for comparison, the one of the confined flare. 
We first notice that the maximum best fit temperature of 
the non-confined flare (~ 4 MK) is much lower than the 
true maximum temperature (see Table ^), by more than a 
factor 2. The reason is that the best-fit temperature is a 
weighted average over the emitting region, which is much 
larger and cooler than the directly heated region; in the 
latter the temperature is close to the plasma maximum 
temperature. We end up with an “observed” temperature 
well below typical flare temperatures, although there is 
plasma above 10 MK and the rate of energy released is 
enough to produce a medium (M GOES class) solar flare, 
as the confined one shown in the figure, with peak plasma 
temperatures around 20 MK. The emission measure val¬ 
ues of the reference model are lower, by a factor 5, than 
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Time (s) 

Fig. 8. Upper panel: Path in the Density-Temperature di¬ 
agram obtained by fitting the ASCA/SIS spectra synthe¬ 
sized from the reference model results with isothermal 
plasma models at various times (s) during the computed 
evolution (solid line). The dashed line has the labeled slope 
and well describes the path of the flare decay. The path 
of the confined flare model is also shown for comparison 
(dotted line). Lower panel: evolution of the temperature 
(solid lines) and of the emission measure (dashed lines) of 
the reference model (thick lines) and of the model of con¬ 
fined flare (thin lines). For the sake of clarity, all quantities 
are normalized to 1: the normalization factors are 4.0 MK 
and 8.9 MK for temperature, respectively, and 0.3 x 10^'^® 
cm“® and 1.4 x 10^® cm“® for emission measure. 

those of the confined model flare, with a maximum value 
of ~ 3 X 10'’’^ cm“® versus ^ 1.6 x 10'^® cm“®. 

The flare path in the diagram is different from the 
one of the confined flare in the heating phase: the 
emission measure reaches very soon its maximum and 


then gradually decreases (leftwards), and the tempera¬ 
ture slightly fluctuates, whereas the model of confined flare 
shows increasing emission measure at constant tempera¬ 
ture (Jakimiec at al. 1992), as also observed in many so¬ 
lar flares (Sylwester et al. 1993). The path is much more 
’’standard” in the decay phase: temperature and emission 
measure both decrease, along the line with positive slope, 
very similar to many observed cases. The value of the 
slope of the non-confined case would suggest a moderately 
heated decay (e.g. Reale et al. 1997, Favata et al. 2000) 
if interpreted as a confined flare. At the end of the decay 
the average temperature stays relatively long at ~ 2 MK. 

The unusual trend of the n-T path in the heating phase 
is directly connected to the relative timing of the evolu¬ 
tion of the density and the temperature (shown in the 
lower panel of Fig. ||). In confined flares, the temperature 
reaches its maximum and begins to decay well earlier than 
the emission measure; the time delay is ~ 200 s in the 
case shown. This does not occur in the reference model 
of non-confined flare (as well as in all the other models 
in Table the emission measure is at its maximum very 
soon, simultaneously to the temperature, and its evolution 
is similar to the temperature evolution and on similar time 
scales. 

In addition to the different timing of the temperature 
and emission measure evolution, also the distributions of 
the emission measure versus temperature during the evo¬ 
lution of the non-confined model flare shows remarkable 
differences from corresponding ones during the evolution 
of the confined model flare. Fig. ^ shows some distribu¬ 
tions of the emission measure versus temperature (EM(T)) 
sampled at various times during the reference model of 
non-confined flare as compared to one EM(T) obtained 
from the confined flare model at the end of the heating 
phase (t = 150 s). In order to evaluate the effect of the flare 
heating on the plasma, the figure shows the EM(T) pro¬ 
files of the reference model subtracted of the unperturbed 
pre-flare EM(T) profile. In agreement to Eig. ||, the max¬ 
imum amount of emission measure is present during the 
heating phase (see profiles at t = 30 s and t = 150 s). The 
effect of the heating is to add emission measure at tem¬ 
peratures above 2 MK, up to ~ 10 MK. Interestingly, the 
heating produces additional, and even larger, components 
of emission measure at lower temperatures. The overall 
EM(T) profile in the heating phase of the reference model 
therefore decreases with temperature, according, approx¬ 
imately, to a power law of index ^ —1.5. This trend is 
qualitatively different from the typical EM(T) profiles of 
flaring (and non-flaring) loops, which instead typically in¬ 
crease with temperature with an index ~ 1.5. As the heat¬ 
ing is switched off, the hottest part of the EM(T) suddenly 
disappears (t = 154 s), and at 200 s there are no pre-flare 
components hotter than 4 MK. Eurthermore, gaps of emis¬ 
sion measure form for T < 1 MK: the plasma in the heated 
region cools and only the plasma in the outward expand- 
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Fig. 9. Distributions of emission measure vs temperature of the reference model at various times during the flare 
evolution: in the heating phase 30 s (histogram), 150 s (solid line)\ in the decay phase 154 s (dotted line), 200 s (dashed 
line). The distributions shown are subtracted by the initial unperturbed distribution at time (t = Os), also shown 
for comparison (dashed-dotted thick line). A distribution obtained from the model of confined flare (at the end of the 
heating phase, t = 150 s) is also shown for comparison (thick solid line). 



Fig. 10. Light curves in the two selected X-ray lines for the reference model and for the confined flare model. 


ing shell contributes to the emission measure in excess of 
the background. 

As shown in Fig. emission measure components are 
present in a wide range of temperatures during the heating 
phase and, therefore, we can expect significant emission in 
the selected spectral lines. Fig. shows the light curves 
of the reference non-confined flare in the selected Mg XI 
at 9.169 A and Fe XXI at 128.752 A lines, assumed al¬ 
ways in equilibrium of ionization. The corresponding light 
curves for the confined model flare are also shown for com¬ 
parison. The light curves of the reference flare globally re¬ 
semble the other synthesized ones at the ASCA/SIS focal 
plane (Fig. Their decay starts simultaneously. While 


the light curve in the hotter Fe XXI line drops by more 
than two orders of magnitude in less than 5 s, due to the 
sudden cooling of the heated region, the Mg XI line decays 
more slowly, because emission measure survives longer at 
temperatures below 10 MK. Once again, the light curves 
look quite different from the corresponding ones of the 
confined flare model. In particular, the cooler Mg XI line 
peaks later than the hotter Fe XXI line, which is more di¬ 
rectly sensitive to the evolution of the heating. The shapes 
of both the line decays for the reference flare are different 
from those of the confined flare, but it may not be trivial to 
discriminate between them in real data, typically coarsely 
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Fig. 11. Fe XXI 128.75 A line profiles synthesized from the reference model and averaged over the heating phase 
(from 0 to 150 s, solid histogram), and over the next 150 s in the decay phase (dashed histogram). Analogous profiles 
from the confined model, averaged over the heating phase (from 0 to 150 s, solid line) and over the next 150 s (dashed 
line) are also shown for comparison. The normalization factors are 3.9 x 10“^®, 6.0 x 10“^^, 8.0 x 10“^^, 2.9 x 10“^® 
erg cm“^ s“^ A“^, respectively. The spectral bin is 0.02 A. The nominal line wavelength is marked (vertical line). 


resolved in time, and sometimes with small signal to noise 
ratio. 


3 . 4 . The other simulations 


As discussed in section 3.2, the evolution of the non- 


As discussed in Sec. 3.1, the basic free parameters of the 


confined plasma has showed the presence of a steady up¬ 
ward moving plasma front, with speed of several hun¬ 
dreds km/s. If favorably oriented with respect to the line 
of sight, we may expect such high speed plasma compo¬ 
nent to produce significant Doppler blue-shifted compo¬ 
nents in resolved lines. Fig. |T^ shows two detailed profiles 
of the Fe XXI 128.75 A line synthesized from the ref¬ 
erence model, in the assumption that the flare is ’’disk- 
centered” with respect to the observer. One profile is av¬ 
eraged over the heating phase (from 0 to 150 s), and the 
other over the next 150 s in the decay phase. The for¬ 
mer profile is practically entirely blue-shifted, while the 
other also shows a significant static component. The blue- 
shifted components both peak at 128.63 A, correspond¬ 
ing to an outward speed component of ~ 300 km/s, and 
are considerably broadened, with a width of 0.05 - 0.10 
A (100 - 200 km/s). Such persistent blue-shifted com¬ 
ponents are not present in the profiles from the model of 
confined flare. The profile averaged over the whole heat¬ 
ing phase is slightly blue-shifted (0.02 A, 50 km/s), due to 
the plasma evaporation, the other profile is practically un¬ 
shifted. Both profiles are much sharper than those of the 
reference model. Indeed considerable blue-shifted compo¬ 
nents (at speed > 100 km/s) are predicted also by models 
of confined flares, but only in the very initial phase of the 
plasma evaporation, generally within the first 60 s in an 
active region loop. Then the blue-shift is washed out by 
the next more gradual, less dynamic (but much brighter) 
evaporation. 


simulations are the pressure at the base of the corona and 
the heating parameters. The explored values are shown 
in Table |^. We have already described in detail the re¬ 
sults of the reference simulation. The plasma evolution 
in all the other simulations invariably resembles the one 
obtained in the reference case, i.e. initial temperature in¬ 
crease and expanding conduction front, subsequent bow 
density front expanding upwards from the chromosphere, 
decay with very fast temperature decrease and continu¬ 
ous expansion of the density bow front. However several 
quantitative differences occur. First of all, the maximum 
temperature is sensitive to the heating parameters and, 
in particular, it increases with the rate per unit volume 
Hq and with the height Zh- Given the high efficiency of 
thermal conduction, the temperature in a flare reaches an 
equilibrium value very rapidly. In steady-state closed coro¬ 
nal structures (loops) the plasma maximum temperature 
is linked to the heating amount by coronal scaling laws 
(Rosner et al. 1978, hereafter RTV). One may then use 
the temperature of a loop to estimate the heating rate 
with the aid of the scaling laws (e.g. Reale et al. 2000, 
Maggio et al. 2000). In the presence of the same physical 
effects at work, we have applied the scaling laws also to 
our configuration despite the scaling laws were devised for 
closed one-dimensional systems. 

By combining the maximum temperature scaling law 
to the volume heating scaling law we obtain: 

where L was the loop half length in the RTV scaling laws, 
T is the plasma maximum temperature and H is heating 
volume rate. If we assume T = T^ax and H = Hq, we 
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invariably obtain L ~ 10® cm for the cases named Ref, 
H30 and HlOO in Table ^ and L ^ 2 x 10® cm for both 
the cases X5H5 and X5H10. From Table |^, we realize that 
these lengths correspond essentially to the “diameter” of 
the heating distributions, i.e. twice the width an- This 
means that in this non-confined configuration the thermal 
distribution of the flaring region is virtually dominated by 
the heating distribution. 


The maximum temperatures of the simulations range 
between 13 and 25 MK. Notice that, in order to produce 
a 25 MK flare, not indeed a major one on the standard of 
stellar flares, we have needed quite a high heating rate of 
2 X 10®® erg/s. We notice also a weak inverse dependence 
of the maximum temperature on the coronal initial base 
pressure: the heating is less effective in an atmosphere 
initially at higher pressure. 


Fig.|^ includes light curves at the ASCA/SIS focal 
plane for all simulations in Table |^, all computed as de¬ 
scribed in Sect. 3.3. Looking at the hgure, the most strik¬ 
ing feature is certainly that, apart from a scale factor 
clearly depending on the amount of energy released in 
the atmosphere, all light curves share a common shape: a 
very fast rise, a relatively constant phase, lasting until the 
heating is switched off, a sudden fast decay by a factor 
2 to 5 on time scale of few seconds, followed by a more 
gradual decay on time scales of few hundreds of seconds. 
This evolution and its timing traces closely the evolution 
of the heating. 

It is straightforward to attribute the fast initial cooling 
to thermal conduction. From EqJ^, the conductive cooling 
time can be approximated as: 


. ngll 


t: 


5/2 


( 11 ) 


where ng is the density in units of 10® cm“^, /g is a char¬ 
acteristic length (10® cm) and Tj is the peak tempera¬ 
ture (10^ K). In practice, the relatively small density value 
makes the cooling by conduction particularly efficient. As 
shown by Table the low density is common to all sim¬ 
ulations and is an intrinsic, important and characterizing 
feature of simulating non-confined flaring plasma, as fur¬ 
ther discussed below. Notice also that the highest pressure 
case (PI) shows a relatively higher density and, in fact, a 
relatively more gradual initial decay. 


Later on, as temperature decreases, cooling by con¬ 
duction becomes less and less efficient and other cooling 
mechanisms come into play, determining the flattening of 
the light curves in the gradual decay phase. In particular, 
in the analysis of the simulation results (Fi gH) we have 
noticed that in the gradual decay phase the flare emission 
originates mostly from the expanding density shell. The 
front propagates in a medium reasonably approximable as 
isothermal, and, therefore, the gradual decay must be at¬ 
tributed to the progressive reduction of the emission mea¬ 
sure. Since, as we have checked, the amount of plasma in¬ 
volved in the density front is approximately (within ~ 20 


%) unchanged with time as the expansion progresses, the 
emission measure in this phase is roughly proportional 
to density. As the front propagates upwards, it maintains 
approximately the same density contrast, as it occurs in 
strong shocks. Therefore the density of the front decreases 
simply because it moves upper and upper through less and 
less dense layers of the stratified corona. The density de¬ 
creases exponentially: 

n ~ no exp(z/sp) no exp(—f/ry) (12) 

where Sp « 2kTf/{mg) ^ 5000 T/ is the pressure scale 
height of the front at temperature Tf for solar gravity, 
and 

Tf = Sp/Vf « 300^ (13) 

is an e-folding time obtained using the fact that the front 
moves at approximately constant speed (Vf ~ 500 km/s). 
The time r/ ~ 300 s, obtained for T/ ~ 3 MK is indeed 
of the order of the decay times reported in last column of 
Table |. 

Of course, the decaying trend changes again as soon as 
any of the above approximations no longer holds, e.g. as 
soon as the shock weakens (as in case PI for t ^ 400 s) or 
is very distant from the origin (as in the case named Long 
for t ^ 900 s). 

All light curves share the same fundamental differences 
from the light curve of the confined model flare: in partic¬ 
ular the decay starts much earlier and starts fast and then 
is more gradual, whereas in the confined flare it starts slow 
and then steepens. It should be noted that, while the only 
way to have a significantly slower flare decay in the non- 
confined atmosphere is to have a heating released over a 
much larger area (Eq. 0), in contrast with flare heating 
theory (e.g. Golub & Pasachoff 1997), the decay of con¬ 
fined flare can be made slower simply considering a longer 
loop, according to the scaling law (Serio et al. 1991): 

wy « 120^ (14) 

4. Discussion 

Modeling a non-confined coronal flare has produced inter¬ 
esting results from the points of view both of pure model¬ 
ing and of diagnostics, especially when compared to con¬ 
fined events. 

We have considered a reference case of a flare triggered 
in an open structure by a heating pulse which would have 
produced a 20 MK flare in a confined structure, namely a 
40,000 km long solar loop. Then we have explored sev¬ 
eral other cases, by increasing the heating rate, there¬ 
fore approching more typical stellar flare conditions, by 
shifting upwards the heating location, as it may occur in 
spatially large stellar coronae, and by considering initial 
atmosphere at higher pressure, which may be present in 
active stars. 

Table ^ shows the main results of the present study, 
by comparing properties of non-confined flare models with 
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Table 2. Relevant properties of flare models® 


Item 

Property 

Confined 

Non-confined 

Diagnostic power: 

tools^ 

Hydrodynamics 

1 

Flaring region 

All loop 

Localized: heated region 

- 

2 

Scale length in RTV seal, 
laws 

Loop half-length 

Diameter of heated region 


3 

Evaporation evolution 

Plasma accumulation 

Expanding fronts 

- 

4 

Decay time scale 

1201/9/\/T7 

early: 

late: SOOVTk.s/Vsoo 

— 

Hydrodynamics & Diagnostics 

5 

Flare temperature 

High= 

Lower 

Low: LRS® 

6 

Plasma density 

High 10“ cm-3) 

Low (;$ 10^° cm“®) 

Low: HRS® 

7 

T-EM relative timing 

Peak delay 

Synchronous 

High: LRS 

8 

Decay 

Gradual 

East-I-Gradual 

Very high: WLC^ 

9 

Line evolution 

Asynchronous 

Synchronous 

?: HRS 

10 

Blue-shift phase 

Short 

Persistent 

Low: Very HRS -f 
orientation® 

11 

EM(T) 

ocr3/2 

oc r-3/2 

?: MRS'* 

12 

Heating tracer 

Hard X-ray l.c. 

Soft X-ray l.c. 

?: WLC 


“ This table summarizes the main differences between the results of confined and non-confined flare models, discussed in more 
detail in Sect. ^ 

^ Level of diagnostic power to discriminate confined and non-confined flare: Low, High, Very high. ”?” means uncertain 
® With the same heating rate: the value depends on the heating rate. 

LRS: Low Resolution Spectra required 
® HRS: High Resolution Spectra required 
^ WLC: Wide-band Light Curve required 

® the flaring region must be located well inside the stellar disk 
^ MRS: Medium Resolution Spectra required 

those of confined flare models and indicating the poten¬ 
tial way of discriminating them through the analysis of 
observations. 

The hydrodynamic evolution obtained from modeling 
non-confined flares is complex and involves many differ¬ 
ent aspects. The model flare is triggered by an impulsive 
heating in localized region of a stratified corona. This heat¬ 
ing produces an immediate but localized temperature in¬ 
crease: proper hot flare conditions occur only within a rel¬ 
atively small region around the heating deposition {Item 
1 in Table |^. This is at variance from a flare occurring 
in confined plasmas, which are typically entirely involved 
in the flare. In our case, thermal conduction has a more 
limited action range: the conduction front rapidly weakens 
due to its expansion and is not able to heat significantly 
plasma very far from the heating deposition place. 

In spite of the dynamic flare evolution, the plasma 
reaches soon a steady state thermal condition, so that the 
maximum temperature and heating rate per unit volume 
are approximately ruled by a scaling law similar to RTV 
scaling laws for coronal loops, in which the heating scale 
place of the loop half-length (Eq. |^, Item 2 

A crucial difference from confined flare events is the 
evolution of the chromospheric evaporation {Item 3 in 
Table ||). In a closed structure, plasma evaporates first 


length takes 
in Table ||). 


rapidly with a strong density front, and then more grad¬ 
ually but continuously until the loop is completely filled 
up and is close to equilibrium conditions. This long-lasting 
evaporation in a confined structure allows density to reach 
high values 10^^ cm“^). The high density (together 
with the flat temperature distribution along the loop) 
makes cooling by conduction relatively inefhcient, lead¬ 
ing to decay times of hours in long loops. What happens 
in the non-confined cases is radically different. There is 
a single evaporation front, which propagates as a strong 
shock continuously outwards from the chromosphere be¬ 
low the heating place, but since plasma is not confined, it 
does not accumulate in spite of further evaporation from 
the chromosphere. In other words, even though the den¬ 
sity in the flaring region soon becomes several times higher 
than it was in the unperturbed condition, it does not grow 
above the strong shock limit (e.g. Rosenau & Frankenthal 
1976), and it never becomes as high as in typical confined 
flares. This is an unavoidable property of a non-confined 
configuration: plasma cannot pile up and therefore can¬ 
not become very dense. This has several implications, but 
probably the most important one, especially from the di¬ 
agnostical point of view, is that, as soon as the heating 
stops, the plasma in the heated region cools very rapidly, 
because thermal conduction is very efficient in a relatively 
low density medium. We have provided a characteristic 








F. Reale et al.: Modeling Non-Confined Coronal Flares: Dynamics and X-Ray Diagnostics 


19 


time scale for this early and fast decay phase (Eq. |ll|) , 
which, for typical conditions, is of the order of few sec¬ 
onds {Item 4 in TaWe[^). 

We have seen that, after the initial drop, the decay 
progresses more gradually, driven by the decrease of den¬ 
sity of the emitting plasma contained within the expand¬ 
ing evaporation front, with a time scale given by the time 
taken by the front to move across the pressure scale height 
(Eq. 1^. However, even the characteristic decay time of 
this phase, invariably not longer than few minutes, is quite 
small as compared to the decay times of typical observed 
large stellar flares (hours or even days), ruled by the con¬ 
duction/radiation cooling time (Serio et al. 1991). 

The thermal and hydrodynamic evolution of the flaring 
plasma leads to a series of diagnostical implications, which 
can be searched for in observations of stellar flares. In the 
first place, we note that non-confined flare require much 
more energy to enhance the plasma temperature {Item 5 
in Table |^) : with the same heating rate both the maximum 
(see Table ^ and the average (see Fig. ||) temperatures of 
the non-conhned flare are nearly halved with respect to 
the corresponding ones of the conhned model. As obvi¬ 
ous consequence one may expect long-lasting stellar flares 
to be also relatively cool, which is opposite to what is 
typically observed, with temperature even reaching 10® 
K (e.g. Maggio et al. 2000). However, since one cannot 
exclude very high energy rates, the high temperatures ob¬ 
served can hardly be taken as a distinctive signature of 
confinement, in the light of our hndings. 

As discussed before, the plasma density should be in¬ 
stead invariably lower in non-conhned flares {Item 6 in 
Table H). This may be directly diagnosed with density- 
sensitive lines detected with high resolution spectrometers 
such as the gratings on board Chandra or XMM-Newton. 
However, it is not trivial to obtain such data for hares and 
to date no such diagnostics is available yet. 

It is instead well-known that in many hares, both so¬ 
lar and stellar, the temperature peaks and starts to decay 
well before the emission measure. This is consistently de¬ 
scribed by models of conhned hares: as the heating drops, 
the temperature begins to decay, while plasma evaporates 
and makes emission measure increase still for some time. 
In our simulations of non-conhned hares, instead tem¬ 
perature and emission measure both decay immediately 
and simultaneously after the heating stops. This diher- 
ence may indeed be a strong indication of conhnement for 
long-lasting stellar hares {Item 7 in Table ^ . 

Another signihcant indication comes from the analy¬ 
sis of the light curves in a wide X-ray passband, such as 
ASC A/SIS {Item 8 in Table In the case of non-conhned 
hares the light curve invariably drops as soon as the heat¬ 
ing stops, with a decay time in the early phase of few 
seconds (Eq. 0). and then to decay more gradually, with 
a time scale of few minutes (Eq. . The early decay of the 
light curve of the conhned hare is much more gradual, and 
becomes even more gradual (hours) if the conhning loop 


is longer than a solar active region loop (Eq. E)- Unless 
we hypothesize unlikely and ad hoc heating depositions, 
e.g. decreasing very gradually and completely driving the 
decay phase, a sudden drop of the X-ray light curve is an 
expected intrinsic and unavoidable feature of flare events 
in a non-confined corona. 

The light curve in selected X-ray lines are also different 
when comparing conhned and non-conhned events {Item 
9 in Table |^). In particular, in the model of conhned hare 
the light curves in two lines at diherent formation temper¬ 
ature are considerably asynchronous: the hotter line peaks 
and decays earlier. In the non-conhned case the two lines 
evolve practically with the same timing. This diagnostic 
may become interesting if hare light curves in single X-ray 
lines are available. 

The presence of long-enduring, high speed, expanding 
plasma fronts may open the opportunity to detect persis¬ 
tent blue-shifted components in high resolution line pro- 
hles {Item 10 in Table H). Uphows of about 400 km/s have 
been seen in the Ca XIX and Fe XXV lines in the early 
thermal phase of some solar hares (e.g. Antonucci et al. 
1987, 1990). Hydrodynamic models of conhned hares pre¬ 
dict that the initial chromospheric evaporation may pro¬ 
duce blue-shifts in X-ray lines, but for very short times 
(less than one minute in active region loops, see Fig. E- 
We have shown that a steady broadened component blue- 
shifted by ^ 300 km/s may be expected in the FeXXI 
line at 128.75 A for the reference non-conhned model. 
This kind of prediction may not provide signihcant indi¬ 
cations when looking at real hare observations for various 
reasons: blue-shifts are expected only if the haring region 
is favourably oriented to the observer; shifts of few hun¬ 
dreds km/s are hardly detectable even with XMM and 
Chandra high resolution grating spectrometers; although 
less persistent, signihcant blue-shifted components may be 
present even during conhned hares in long loops. All these 
limitations may make any diagnostics of conhnement from 
line shifts rather difhcult. 

Recent works have pointed the attention to the deriva¬ 
tion of the distribution of emission measure with temper¬ 
ature {Item 11 in Table ^ from the analysis of X-ray stel¬ 
lar data (e.g. Drake et al. 2000, Linsky & Gagne 2001). 
Non-conhned hare models predict a hare emission mea¬ 
sure distribution monotonically decreasing with tempera¬ 
ture during the heating phase and then losing the highest 
temperature components and some intermediate temper¬ 
ature components in the decay phase. These distributions 
are quite different from those typically derived from har¬ 
ing and non-haring solar (e.g. Peres et al. 2000, Reale et 
al. 2001) and non-haring stellar data, which instead show 
well-dehned EM peaks. This may be a discriminating fea¬ 
ture, whenever medium resolution spectra provide reliable 
emission measure distributions during stellar hares (e.g. 
Osten et al. 2000). 

As a hnal diagnostic consideration, we point out that, 
if data show evidence of a hare occurring in a non-conhned 
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plasma, then according to our modeling the soft X-ray 
light curve would faithfully trace the evolution of the heat¬ 
ing, an important diagnostics for coronal physics [Item 12 
in Table An equivalent information is generally derived 
from hard X-rays in solar flares (e.g. Golub & Pasachoff 
1997). 

We may wonder whether our results, such as the light 
curves presented in Fig. are affected by the Non Equi¬ 
librium of Ionization (NEI) conditions of the plasma. In 
general, in fact, the X-ray emissivity depends on the ther¬ 
mal history of the emitting plasma parcel, which is un¬ 
fortunately not straightforward to evaluate. Bocchino et 
al. (1997) show that the 0.1-2.0 keV X-ray emissivity of 
an impulsively heated plasma (a situation similar to the 
rise of the flare in our models) approaches the equilibrium 
value for ionization time logr ^ 11 in s cm“^. This is also 
true in the ASCA/SIS bandwidth. In the initial phase of 
the flare the typical density of the bright X-ray emitting 
plasma is in excess of 10® cm“^, corresponding to an equi¬ 
librium time scale of the order of 100 s: we may expect a 
modification of the light curve before 100 s, consisting in a 
higher NEI value reached at the end of the rising phase and 
a slow descend to the sustained rate. In the fast decay at 
the end of the heating phase the NEI conditions may have 
some importance and bring to variations of the light curve 
in this phase. However, these effects will be reduced for 
hotter and hotter events, i.e. typical intense stellar flares, 
because the continuum emission dominates more and more 
on the line emission. Later on, in the gradual decay, most 
of the emission comes from the expanding shock, in which 
we have only a slight decrease of the temperature, and 
therefore the plasma is only marginally overionized. The 
small deviation from equilibrium gives quasi-equilibrium 
conditions during most of the decay phase, thus we do not 
expect substantial modifications of the global time scales 
we have worked out. 

In summary, there are two basic features that may 
characterize invariably flares in non-confined regions and 
may be relatively easy to trace in stellar X-ray data: 
light curves with very fast decays and synchronous evo¬ 
lution of density and temperature. Both these features 
are in contrast with evidence from stellar flares. Then we 
come to an interesting conclusion: the long-lasting X-ray 
flares typically detected in several stellar observations can¬ 
not occur in non-confined atmospheres. This is somewhat 
unexpected, because long-lasting stellar flares have been 
typically associated to very large coronal structures (see 
Sect. 1^, which have the non-confined configuration as an 
asymptotic extreme. This work shows that most likely the 
long-lasting flares, such as those in active stars, still oc¬ 
cur in closed, although large, structures and that the role 
of confinement of the coronal magnetic field must be in¬ 
variably significant in those events. In other words, the 
magnetic fields must always be strong and/or the mag¬ 
netic structure never breaks open. In the same direction, 
notice also that non-confined flares appear to involve huge 


amounts of heating (see Table ||) which may not be en¬ 
tirely realistic when compared to the energy budget of 
stellar coronae. 

If long-duration stellar flares are unlikely to be de¬ 
scribed as breaking the magnetic confinement, are there 
any other observed events which instead could be? As 
mentioned in Sect. long stellar flares are preferentially 
detected in X-ray observations. However there is a class 
of coronal variations which occur on small time scales 
and which have been detected on dMe star UV Ceti by 
ROSAT (Schmitt et al. 1993). The possibility that such 
short events may be interpreted as small but very intense 
non-confined events should be explored. 

This study in some way shows an interesting theoret¬ 
ical perspective to be compared with SNR models and 
Coronal Mass Ejections (CME) models. Our model shows 
an overheated plasma which expands under the effect of 
heating in a hot and thermally conducting stratified at¬ 
mosphere, and we provide for it characteristic scalings of 
general validity (Eqs. |^, |^). For instance, since the 

time scale of the late decay tj (Eq. |^) depends practi¬ 
cally only on the speed of the evaporation front and on the 
thermal plasma conditions far from the flaring region, it is 
little dependent on the flare heating and on other model 
parameters: it is a general decay time scale for shell fronts 
expanding in a stratified conducting corona. It should be 
noted that application to other astrophysical systems is 
not trivial. In SNR models, plasma expands in a much 
cooler, thinner and less conducting medium. On the other 
hand, CME’s are mostly observed in the UV band, and 
are therefore to be modeled as relatively cool perturba¬ 
tions, whereas the non-confined fronts modeled here are 
at standard coronal temperatures. For the latter fronts to 
become proper CME’s, a mechanism to thermally insulate 
the fronts from the surrounding corona should be invoked, 
so that they are free to cool by radiation and to emit in 
the UV band (e.g. Ciaravella et al. 2001). 

5. Conclusions 

Stellar flares are generally observed to evolve and decay 
on time scales ranging from several hours to days, suggest¬ 
ing large flaring regions or even lack of magnetic confine¬ 
ment. This work shows that the hydrodynamic evolution 
of flares occurring in non-confined atmospheres leads in¬ 
variably to a much faster decay of the brightness (on time 
scales of very few minutes) after the heating phase, in a 
wide range of the physical parameters. One of the main 
implications is that the long duration of stellar flares is 
not indicative of non-confined plasma. Reversing the argu¬ 
ment, it is highly probable that the observed long-lasting 
stellar X-ray flares involve forms of plasma confinement in 
closed coronal structures. We have pointed out other char¬ 
acteristic features, with more or less diagnostical power, 
e.g. the synchronous evolution of density and temperature, 
summarized in Table that seem to indicate that X-ray 
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observed stellar flares involve mostly plasma confined in 
closed structures. 

Although this work simulates rather extreme condi¬ 
tions of flares occurring in completely non-confined atmo¬ 
spheres, the result that long-lasting stellar flares are likely 
not occurring in large open structures largely motivates 
its validity. 

Furthermore this work predicts the possible existence 
of a new phenomenological class of events, characterized 
by peculiar light curves with a limited range of time scales. 
Such kind of events may be addressed by observations at 
high sensitivity levels, which allow for high time resolution 
and high signal-to-noise ratio, such as those obtained from 
the current missions Chandra and XMM-Newton. 
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